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Preface 


The  study  reported  herein  was  conducted  as  part  of  the  Chesapeake 
Bay  Three-Dimensional  Model  Study.  It  was  sponsored  by  the  Chesa¬ 
peake  Bay  Program  Office,  U.S.  Environmental  Protection  Agency,  and 
the  U.S.  Army  Engineer  District,  Baltimore. 

The  principal  author  of  this  workplan  was  Dr.  Carl  F.  Cerco,  Water 
Quality  and  Contaminant  Modeling  Branch  (WQCMB),  Environmental 
Laboratory  (EL),  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES).  The  chapter  on  Hydrodynamic  Model  Activities  was  formulated 
in  cooperation  with  Dr.  Billy  Johnson,  Waterways  Division,  Hydraulics 
Laboratory,  WES.  Technical  review  was  provided  by  Mr.  Ross  Hall, 
WQCMB,  and  Dr.  Mark  Dortch,  Chief,  WQCMB. 

Management  was  provided  by  Mr.  Donald  L.  Robey,  Chief,  Environ¬ 
mental  Processes  and  Effects  Division,  EL.  Overall  supervision  was  pro¬ 
vided  by  Dr.  John  W.  Keeley,  Director,  EL, 

At  the  time  of  publication  of  this  report,  WES  Director  was  Dr.  Robert 
W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 

This  report  should  be  cited  as  follows: 

Cerco,  D.  F.  (1994).  “Workplan  for  tributary 
refinements  to  Chesapeake  Bay  eutrophication  model 
package,”  Miscellaneous  Paper  EL-94-S,  U.S.  Army 
Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 


Chapter  1 
Introduction 


The  Corps  of  Engineers,  in  partnershq)  with  the  USEPA  Chesi9)eake  Bay 
Program  Office,  recently  completed  a  diree-dimensional  modd  study  of 
eutrqriiication  in  Chesiq»ealce  Bay  and  tributaries.  The  model  package  sillied 
included  an  intratidal  hydrodynamic  model  (Idinson  et  al.  1991),  an  intotidal 
water-quality  model  (Cerco  and  Cole  1994),  and  a  bendiic  sediment  diagenesis 
model  (DfToro  and  Fhapatrick  1993).  The  Chesq>eake  Bay  Program  Office 
has  prt^iosed  a  series  of  tributary  rdiimneitts  to  &e  existing  model  package. 
The  refinements  have  die  tbllowing  goals: 

To  provide  a  more  accurate  simulation  of  the  fitte  and  transport  of 
nutrients  in  die  tributaries. 

To  simulate  the  effects  of  light  attenuation  by  suspended  sediment  and 
nutrients  on  submerged  aquatic  vegetation. 

To  more  accuratdy  simulate  dissolved  oxygen  concentrations  in  die 
tributaries. 

To  provide  the  basic  modeling  framework  for  future  toxics  modding. 

To  provide  a  more  detailed  hydrodynamic  structure  diat  can  be  qiplied 
to  the  simulation  of  oil  spills  in  the  James,  Yoric,  and  Riqipahannock 
Rivers  and  Baltimore  Harbor. 

The  mcisting  roodd  provides  excdlent  rqiresentation  of  physical  and 
eutrc^hication  processes  in  the  mainstem  bay.  Rqiresentation  in  the 
tributaries  is  uneven,  however.  A  broad  characterization  is  diat  modd 
p^formance  diminishes  aloi^  widi  die  size  of  die  tributary  and  according  to 
distance  from  the  juncture  widi  die  mainstem.  Consequendy,  representation 
of  the  larger  tributaries,  die  James  and  die  Potomac,  is  siqierior  to  lesser 
tributaries  induding  die  York,  Ranishannock,  and  Patuxent.  Modd 
performance  in  die  lower  Potomac,  adjacent  to  the  Bay,  is  superior  to 
performance  in  die  freshwater  Potomac,  immediately  below  die  fidl  line.  In 
dl  cases,  the  existing  modd  is  unable  to  accuratdy  quantify  water-qudity 


benefits  of  a  pn^wsed  forty^)CTcent  reduction  in  controllable  nutrient  loads. 

Starting  points  for  improved  rqpresmitation  in  the  tributaries  have  been 
identified  (Cerco  and  Cole  1994).  They  are: 

Magnitude  and  composition  of  loads  must  be  quantified  with  extreme 
accuracy. 

nienomena  near  the  limit  of  salt  intrusion  must  be  rq)resented.  These 
include  flow  convergence,  accumulation  of  organic  a^  inorganic 
particles,  and  occurrence  of  a  chlorophyll  maximum. 

Longitudinal  and  lateral  resolution  of  the  tributary  grids  should  be 
improved. 

At  present,  living  resources  are  indirectly  coupled  to  the  model.  Model 
predictions  of  dissolved  oxygen,  nutrient  concentrations,  and  li^t  extinction 
are  compared  to  "living-resources  criteria".  The  comparison  indicates  die 
effect  of  management  activity  on  conditions  favorable  to  living  resources  but 
does  not  quantify  the  living-resource  response.  Nor  does  the  current 
framework  compute  feedback  effects  in  which  living  resources  impact  water 
quality.  At  present,  only  the  effects  of  water  quality  on  living  resources  are 
quantified. 

The  present  document  comprises  a  worlqilan  to  inqirove  tributary 
rqiresentation  and  to  incorporate  living  resources  direedy  into  the  model 
framework.  Four  tributaries  have  been  sdeaed  by  the  Chesqieake  Bay 
Program  Office  for  emphasis  under  this  tributary  refinements  program.  They 
are  the  James,  York,  Rtqipahannock  Rivers,  and  Baltimore  Harbor.  The 
James,  York  and  Rsqipahannock  w^e  specified  because  tributary-specific 
models  are  required  to  address  water-quality  and  living-resource  bmiefits  to  be 
derived  from  nutrient  reduaions.  Baltimore  Harbor  was  specified  because  it 
presents  unique  management  problems,  coupled  with  long-term  toxics  inqiacts, 
which  cannot  be  addressed  in  the  current  model  framework.  The  time  scale 
for  the  project  is  four  years  from  initiation  to  completion.  Anticipated 
commencement  is  April  1,  1994. 


1-2 


CiMparl  knrodiietioii 


Chapter  2 

Woiicshop  Summary 


Three  technical  workshc^  were  oonveiMd  to  assist  in  formulation  of  a 
technical  ^>proadi  for  tributary  refinonents.  The  woricshops  were  attended 
by  memben  of  government,  anemia,  and  private  enterprise  widiin  and 
outside  the  Bay  community.  Conclusions  of  foe  worksh(^  are  summarized 
below.  Complete  proceedings  are  available  in  HydtoQual  (1993). 

Workshop  I:  Tributary  Models  and  Main  Bay  Model 
Linkage 

The  workshop  was  convened  primarily  to  investigate  mefoods  by  wfaidi 
finorscale  tributary  models  will  interface  with  a  larger-scale  model  of  foe 
mainstem  bay.  Among  foe  topics  discussed  were  grid  resolution  in  foe 
tributaries,  interactions  between  foe  tributaries  and  mainstem,  and  target 
processes  to  be  rq)roduced  in  foe  model.  Significam  recommendations 
included: 

Interactions  between  foe  mainstem  and  tributaries  indicate  foat  foe 
tributaries  cannot  be  modelled  indq>endently  of  foe  mainstem. 

The  hydrodynamic  model  should  accurately  rq[>resent  bafoymetry  to 
foe  2  m  contour,  at  least. 

Spatial  resolution  of  noqmint-source  loads  to  foe  water-quality  modd 
should  be  refined. 

Target  processes  to  be  rqiresrated  in  foe  modd  include:  tributary 
turbidity  maximum,  tributary  chlorophyll  maximum,  tilting  evoits  at 
foe  bay-tributary  interfaces,  and  tiir^ayer  circulation  in  Baltimore 
Harbor. 

Modd  plication  should  be  extended  beyond  foe  current  1984-1986 
period. 
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The  formulation  of  the  opmi-mouth  boundary  condition  should  be 
reviewed.  Extension  of  the  grid  outside  foe  bay  moufo  should  be 
considered. 

Workshop  II:  Living-Resource  Processes  and  Ecosystem 
Modeling 

The  workshop  addressed  foe  objectives  in  adding  living  resources  to  foe 
model  framework.  Next,  state  variates  and  processes  for  addition  to  foe 
existing  water-quality  modd  were  suggested.  Significant  recommendations 
included: 

The  hydrodynamic  model  should  resolve  bafoymttry  to  foe  2  m 
contour.  Volume  within  foe  2  m  contour  should  be  modeled  as  a 
littoral  zone  treated  as  storage  area  by  foe  hydrodynamic  model. 

Mechanistic  representation  of  submerged  aquatic  vegetation  should  be 
added  to  foe  model  package.  The  SAV  submodel  should  include  roots 
and  rhizomes,  leaves,  epifauna,  and  grazers  on  epifouna. 

Mechanistic  rq)resentation  of  bathos  should  be  added  to  foe  sediment 
submodel.  Rq)resentation  of  benthos  requires  modeling  of  benfoic 
algae,  dissolved  organic  carbon,  and  two  groups  of  macro¬ 
invertebrates. 

Zooplankton  should  be  added  to  foe  water-quality  model. 

Continuous  salinity  records  should  be  collected  at  midchannel  and  in 
shoal  areas  of  tributaries.  The  records  will  be  used  to  verify  modelled 
exdiange  b^een  foe  littoral  zone  and  cbannd  conveyance  area. 

Workshop  III:  Suspended  Sediments  and  Light 
Attenuation 

The  workshop  considered  foe  addition  of  suspended  sediments  to  foe 
model  framework.  The  potential  of  a  fullyi>redictive  q>proach  and  of 
alternatives  wae  considered.  Effects  of  suspended  sediments  on  light 
attenuation  were  examined  and  inclusion  of  foe  effects  in  foe  model  were 
discussed.  Significant  recommendations  wo’e: 

A  folly-predictive  suspended-sediment  model,  operating  on  foe 
temporal  and  spatial  scales  of  foe  hydrodynamic  model,  is  desirable 
but  not  immediately  feasible. 
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The  first  step  in  sediment  modeling  requires  investigation  of  die 
sedimmit  budget  of  the  bay.  The  existing  model  framework  should  be 
used  for  this  investigation. 

Light  attenuation  can  be  successfully  computed  as  a  function  of 
organic  and  inorganic  solids,  dissolved  carbon,  and  chlorophyll  *a’. 

Successful  modeling  of  SAV  in  the  littoral  zone  requires  consideration 
of  interactions  between  suspended  sedimrats  and  light  attenuation. 
Sediment  resuspension  by  wind  must  be  included  in  the  representation. 

Consideration  should  be  given  to  special  studies  to  determine  the 
relationship  between  light  and  solids  and  to  measure  sediment 
resuspension  as  a  function  of  wind  velocity. 
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Chapter  3 

Hydrodynamic  Model  Activities 


The  Hydrodynamic  Model 

Eii^>loyment  of  die  CH3D-WES  (Qin^utadonal  Hydrodynamics  in  Three 
Dimensions  •  Waterways  Expmment  Station)  hydrodynamic  model  wUl 
continue  during  tributary  refinemaits.  Modd  formulation  is  based  on 
principles  expressed  by  the  equations  of  motion,  conservation  of  volume,  and 
conservation  of  mass.  Quantities  computed  by  dir  modd  include  three* 
dimensiond  vdochies,  surface  devadon,  voticd  viscosity  and  diffosivity, 
temperature,  sdinhy,  and  density.  Details  of  the  modd  formulation  and  prior 
^plication  to  Chesapeake  Bay  are  presented  by  Johnson  et  d.  (1991). 

Model  Modifications 

The  CH3D-WES  code  will  be  modified  to  incorporate  a  littord  zone 
around  the  shordine  of  the  bay  and  tributaries.  Tie  littord  zone  will  extend 
from  the  land-water  interface  to  the  two-meter  (at  mean  tide)  depth  contour 
(Figure  3-1).  Tlie  "littord  zone”  is  a  new  name  for  a  distinction  betweoi 
storage  and  conveyance  area  that  has  been  onployed  for  ova  twenty  years 
(Thatcher  and  Harleman  1972)  and  is  still  us^  in  representing  die  geometry 
of  bay  tributaries  (Park  and  Kuo  1993).  The  storage  area  fills  and  entities  as 
surface  levd  rises  and  fdls  at  a  diannd  cross  section.  Laterd  exdiange 
between  die  storage  and  conveyance  area  is  computed  according  to  continuity. 
Axid  flow  occurs  in  die  conveyance  area  only  a^  is  computed  according  to 
the  coupled  momentum  and  condnuity  equations.  The  momoitum  equation 
employs  die  cross  section  of  die  conveyance  area  while  die  continuity  equation 
enqiloys  die  totd  cross  section  of  die  conveyance  and  storage  areas. 

A  second  modification  to  (^3D-WES  will  dlow  employment  of  only  one 
cdl  in  die  verticd.  At  present,  two  cells  are  required  idiid  restricts  the 
minimum  dqidi  rqiresented  in  die  modd.  Enqiloyment  of  one  verticd  cdl 
dqnh  and  addition  of  die  littord  zone  will  provide  excdlent  rqiresoitation  of 
geometry  and  flows  throughout  die  diannd  cross  section. 
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Computational  Grid 

Virginit  TributiriM 

A  unique  system-wide  grid  will  be  constructed  for  eadi  Virginia  tributary. 
On  each  grid,  present  resolution  widiin  one  tributary  wOl  be  doubled  along  foe 
longitudinal  and  lateral  axes.  Existing  vertical  resolution  (surfoce  layer  «>2m, 
all  others  « l.Sm)  will  be  maintained.  Outside  foe  tributary,  resolution  will 
gradually  transition  to  an  equivalent  of  foe  present  grid.  The  strata  offers 
two  advantages: 

Resolution  in  foe  tributaries  is  enhanced  vfoile  total  conq>utational 
burden  is  minimized. 

Interactions  of  foe  tributary  wifo  foe  mainstem  and  ofoer  tributaries 
are  rq)resented. 

Present  grid  size  and  tqtproximate  proposed  grids  are  summarized  in  Table 
3-1. 


Table  3-1 

Existing  and  Proposed  Grid  Size 

ExMnf 

Smlice  Ceils 

EidstiRi 

Tooil  Ceils 

Rroposei 
Smfece  Cells 

Propesei 

Telal  Cells 

James  River 

42 

139 

168 

556 

York  River 

61 

224 

244 

896 

Rappahannock 

River 

35 

120 

140 

480 

1  System  Total 

729 

4,029 

•  1300 

•  9,000 

1 

1  Baltimore 

I  Harbor 

15 

76 

150 

780 

Baltimore  Harbor 

Baltimore  Harbor  (including  foe  Patapsco  and  Back  Rivers)  is  less 
extensive  foan  foe  ofoer  tributaries  select^  for  refinements.  Detailed 
resolution  of  geometry  in  foe  harbor  means  foat  grid  cells  will  be  much 
smaller  foan  elsewhere  in  foe  system.  The  diminutive  cells  presoit  two 
problems.  The  first  is  transition  between  foe  harbor  and  adjoining  bay. 

Large,  discontinuous  junq>s  in  grid  size  are  not  permissible.  Consequently, 
foe  fine  grid  will  extend  into  foe  bay  and  infiuence  grid  scale  forou^ut  foe 
system.  The  second  problem  is  numerical  time  step.  The  time  stq>  forced  by 
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the  fine  grid  spacing  in  the  harbor  will  determine  die  time  stqi  throughout  die 
syston.  The  solution  is  to  create  a  fine-grid  modd  of  die  haiiior  and  adjacent 
bay  which  will  operate  indqiendendy  of  the  system-wide  modd.  Approximate 
extent  of  die  grid  will  be  from  the  Susqudianna  fall  line  to  die  Bay  Bridge, 
near  Amugmlis.  Information  from  hyd^ynamics  runs  on  one  of  die  system- 
wide  grids  will  be  used  to  set  boundary  conditions  at  the  (^len  end  of  the 
Baltimore  harbor  grid. 

Remainder  of  the  System 

Regridding  the  Virginia  tributaries  will  force  grid  refinmnents  throughout 
the  system.  System-wide  refinements  are  required  to  avoid  large 
discontinuous  jumps  in  grid  scale  from  die  tributaries  to  the  bay.  Also, 
gridlines  must  be  continuous  from  one  boundary  to  the  other.  Consequendy, 
lateral  grid  resolution  at  die  mouth  of  die  Rappahannock,  for  example,  must 
be  carried  across  the  mainstem  to  the  Eastern  Shore. 

As  part  of  the  regridding  process,  refinements  will  be  made  in  die 
Potomac  to  rectify  some  existing  problem  areas.  Notably,  die  right-angle 
bends  at  Mathias  Point  and  Maryland  Point  will  be  smoofoed.  The  littoral 
zone  wUl  be  installed  throughout  die  system  including  die  Potomac, 
Susqudianna  Flats,  and  Tangier  Sound. 

Open-Mouth  Boundary  Conditions 

Good  practice  in  predictive  modeling  is  to  place  boundary  conditions  as 
far  as  possible  from  regions  of  interest.  Downstream  boundaries  in  the 
current  grid,  at  the  mouth  of  die  bay,  are  close  to  die  James  River,  a  primary 
focus  of  the  tributary  refinements.  Prospects  for  extending  the  boundary  away 
from  the  James,  out  into  the  continental  shelf,  are  cloudy,  howevw.  Among 
the  concmns  are: 

Availability  of  hydrographic  and  water-quality  observations. 

Computational  problems  associated  with  specification  of  three  open 

boundaries,  and 

Additional  computational  burden. 

Opinions  at  the  workshops  were  mixed  as  to  die  feasibility  and  desirability 
of  extending  die  open  boundary.  At  initiation  of  the  study,  a  careful 
examination  of  die  above-noted  issues  will  be  conducted  before  determination 
of  the  downstream  boundary  location.  Examination  will  include  an 
exploratory  hydrodynamic  model  run  with  die  grid  extended  *30  km  outside 
the  bay  mouth.  Dmation  of  die  run  will  be  thirty  to  sixty  days.  One 
objective  is  to  determine  whedier  observed  conditions  at  the  bay  mouth  can  be 
rqiroduced  when  boundary  conditions  are  specified  on  the  exteiided  grid. 
Model  performance  at  the  bay  moudi  will  be  a  primary  considoation  in 
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determiiiatkni  of  die  locttkm  of  die  grid  boundary. 

Stqutnee  to  be  Modeled 

Vlr^aia  Tritataiies 

Hie  York  River  is  ^ecified  as  the  first  to  be  modeled  since  h  contains 
extensive  littoral  r^ions  and  living  resources.  Model  runs  tm  the  Y<xk  grid 
will  provide  die  tasting  ground  for  code  revisions  and  implenicnration  of  foe 
litttxal  xone.  The  James  and  Rj^ipahannodc  will  be  modeled  simultaneously 
fUlowing  initial  efforts  in  foe  Yoric. 

Periods  sdected  for  modd  calibration  widiin  the  tributaries  will  be 
qiecified  following  examination  of  availdile  data.  After  calibration, 
production  runs  will  be  conducted  for  1985-1987  and  1994.  Additional 
verification  of  foe  modd  will  be  conducted  usipg  sdiniQr  data  fttmi  the 
production  years.  Hydrodynamic  production  runs  cmre^ond  to  years  sdected 
for  cdibration  and  verification  of  foe  water-qudity  modd. 

BaMmoro  Harbor 

Woric  on  foe  indqiendent,  Bdtimore  Harbor  modd  wfll  commence 
following  cdibration  of  foe  Virginia  tribittaries.  Ddaying  Bdtimore  Harbor 
until  Year  2  of  foe  project  will  dlow  for  utilization  of  a  system-wide  grid  to 
set  downstream  boimdary  conditions  for  the  fine-grid  harbor  modd. 

Periods  sdected  for  cdibration  of  foe  harbor  modd  will  be  qwcified 
following  examination  of  avaflable  data.  Timdy  cdibratkm  and  initid 
verification  of  foe  modd  require  enqiloyment  of  data  collected  prkxr  to  die 
1994-1995  surveys.  After  (Vibration,  production  tuns  will  be  conducted  for 
1985-1987  and  1994-1995.  Additiond  verification  of  foe  modd  will  be 
conducted  using  sdinity  data  from  foe  production  years.  Hydrodynamic 
production  runs  correHiond  to  years  sdected  for  cdibration  and  verificatkm  of 
foe  water-qudity  modd. 

We  note  that  cdibration  to  1$194-1995  within  foe  time  frame  of  the  project 
will  be  problematic  and  may  be  inqiossible.  Cditoation  requires  obser^ 
sdinity  and  tidd  boundary  conditions,  meteorologic  conditions,  and  freshwater 
runoff  records.  These  mi^  not  be  available  wifoin  foe  time  fri^  of  die 
project.  If  foe  system-wide  modd  is  to  be  used  to  provide  boundary 
conditions,  an  additkmd  year,  1995,  nust  be  added  to  foe  four  yean  of 
planned  production  runs.  Data  for  foe  system-wide  runs  may  not  be  availaUe 
dfoer. 

System-Wide  Model 

A  single  system-wide  grid,  incorporating  dl  refinements  to  Virginia 


tributaries,  may  be  created  following  hydrodynamic  model  calibration  on 
individual  grids.  Production  tuns  on  this  grid,  198S>1987  and  1994,  will  be 
completed  if  required  to  accurately  represent  tributary  interactions  (hiring 
nutrient  load-reduction  scenarios.  Determination  of  die  need  for  a  single  grid 
will  be  made  following  transport  tests  on  individual  grids. 

Data  Bases 

Virginia  Tributariea 

Extensive  observations  of  tides,  currents,  and  salinity  have  been  collected 
in  the  Virginia  tributaries  by  the  Virginia  Insthute  of  Marine  Science  (VIMS). 
Cooperation  widi  VIMS  will  be  solicited  in  enqiloying  diis  data  base  for 
model  calibration.  We  wUl  also  seek  advice  for  target  processes  diat  should 
be  simulated,  such  as  spring-neap  stratification  cycling  and  tributary-mainstem 
exchange  processes.  The  VIMS  data  will  be  supplemented  widi  salinity 
observations  collected  in  the  Chesapeake  Bay  Program  monitoring  program. 

Bahimere  Harbor 

Our  workshops  revealed  limited  hydrogrqihic  data  exists  for  Baltimore 
Harbor.  The  only  known  current  (foservations  were  collected  in  1979. 
Monitoring  program  stations  are  limited  to  one  eadi  in  the  Pat^Mco  and  Back 
Riv<»s.  These  data  will  be  enq)loyed  in  initial  development  of  the 
hydrodynamic  model.  Hydtogrtqihic  data  collet  by  die  Maryland 
Elqiartment  of  the  Environment  in  1994  will  be  incorporated  into  the 
hydrodynamic  model  data  base  as  it  becomes  available. 

Linkage  to  the  Water  Quality  Model 

In  the  Chesiyieake  Bay  Three-Dimensional  Model  Study,  Lagrangian 
averaging  (Dorttdi  1990)  was  permed  on  hydrodynamics  prior  to 
employmmit  in  die  water  quality  modd.  A  processor  imbeclded  in  the  CH3D- 
WK  code  transformed  vdocities,  surface  devations,  and  verticd  diffosivities 
computed  on  a  five-minute  basis  into  tidd-average  vdues  ouqmt  every  12.5 
hours.  Employment  of  die  Lagrangian  dgoridim  reduced  die  storage 
requirements  for  the  hydrodynamic  ouqiut  by  an  order  of  magnitude.  Expert 
consensus  is  the  Lagrangian  dgoridim  is  not  suited  for  use  in  die 
geometricdly-confined  tributaries.  Consequendy,  intratidd  (<  12.4  hours) 
radier  than  intertidd  (>  12.4  hours)  hydrodynamics  will  be  used  to  drive  foe 
water  qudity  model. 

Averaging  hydrodynamics  up  to  a  time  scde  more  lengfoy  dian  foe 
CH3D-WES  time  stqi  should  still  be  possible,  however.  A  target  of  one-  to 
two-hour  hydrodynamic  increments  is  desirable.  Intratidd  averaging  of 
hydrodynamic  ouqiut  is  enqiloyed  to  minimize  foe  storage  requirement  of 
hydrodynamic  iiqmt  to  foe  water-qudity  modd.  The  water-^idity  modd  time 


step  is  not  necessarily  e^ial  to  the  hydrodynamic  averaging  increment.  The 
typical  water-quality  model  time  stq)  enqUoying  the  Lagrugian-average 
vdocities  was  two  hours.  Larger,  intratidal  vdocities  nd  nndler  grid 
spacing  will  likdy  force  wder-qudhy  modd  time  steps  mudi  less  dian  two 
hours  during  tributary  refinonents. 

During  the  Three-Dimensional  Modd  Study,  extensive  experience  was 
gained  in  coupling  hydrodynamic  and  watw-quality  modds  and  in  testing  die 
linkage.  Tests  focus^  on  comparisons  of  transport  of  ctmservative  substances 
in  both  modds.  Initial  tests  oonqiated  "dye  dunqis*  in  the  hydrodynamic  and 
water  quality  modds.  Secondary  tests  were  examinations  of  sdinity  conqnited 
in  both  modds.  SimUar  testing  will  be  conducted  during  tributary  refinemems 
to  oisure  correct  and  accurate  transfer  of  information  from  die  hydrodynamic 
to  water-quality  modds. 
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Chapter  4 

Water-Quality  Model  Activities 


The  Water-Quality  Model 

Enqiloyiiiieiit  of  the  CEK^UAL-ICM  water-tpiality  modd  will  coadmie 
during  tributary  refinements.  Tbe  fbundatkm  of  CE-QUAL-IO^  is  the 
solution  to  die  three-dimensional  mass-conservation  equation  for  a  ocmtrol 
volume.  Control  volumes  in  CE-QUAL-ICM  correqiond  to  cdls  in  x-y-z 
qiace  (m  die  CH3D-WES  grid.  CE-QUAL-ICM  solves  die  conservation  of 
mass  equation  in  eadi  volume  for  eadi  modd  state  variaUe. 

A  suite  of  twenty-two  state  variables  (Table  4-1)  is  currendy  employed  to 
modd  eutrcqdiication  processes  in  die  water  cdumn.  Kinetic  interactions 
afiecting  die  state  variables  are  described  in  over  80  paitid  difforentid 
equatioiiis  diat  require  evahu^n  of  over  140  paraiiie^.  Ibe  kinetics 
describe  carbon,  phoqihorus,  nhrogeii  and  silica  cydes  and  die  dissolved 
oxygen  bdance.  Deta^  of  modd  formulation  and  prior  ^iplication  to 
Chesqieake  Bay  are  found  in  Cerco  and  Cole  (1994). 

The  Sediment  Model 

A  unique  feature  of  the  eutroidiication  modd  package  is  coupling  of  die 
water-quality  modd  with  a  predictive  modd  of  diagenesis  in  benthic 
sediments.  The  modd  is  driven  by  net  settling  of  organic  matter  from  die 
water  column  to  the  sediments.  In  the  sediments,  die  modd  simulates  die 
diagenesis  (decay)  of  die  organic  matter.  Diagenesis  produces  oxygen  demand 
and  inorganic  nutrients.  Oxygen  demand,  as  sulfide  fin  sdtwater)  or  meriiane 
On  freshwater),  takes  diree  padu  out  of  die  sediments;  export  to  the  water 
column  as  di^cd  oxygen  demand,  oxidation  at  die  sediment-water  interfile 
as  sediment  oxygen  demand,  ot  burid  to  deep,  inactive  sediments.  Incvganic 
nutrients  produced  by  diagenesis  take  two  paths  out  of  the  sediments:  release 
to  die  water  column,  or  burial  to  deep,  inactive  sediments.  A  listing  of 
sediment  modd  state  variables  and  p^icted  sediment-water  fiuxes  is  povided 
in  Table  4-2.  Detafls  of  modd  formulation  are  found  in  DiToro  and 
Fitzpatrick  (1993). 
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TabI*  4-1 

Water  Quality  Modal  State  Variablaa 

Tomporatura 

Saiinihf 

Total  Activo  Motal 

Cyanobacteria 

Diatoms 

Green  Algae 

Oissolwd  Organic  Carbon 

table  Particulate  Organic  Carbon 

Rafractory  Partkulata  Organic  Carbon 

Ammonium 

Nitiata 

Dissolved  Organic  Niuogen 

Labiia  Partkulati  Organic  Nitrogen 

Refractory  Particulate  Organic  Nitrogen 

Total  Phosphate 

Oissofred  Organic  Phosphorus 

Labile  Particulate  Organic  Phosphorus 

Refractory  Particulate  Organic  Phosphorus 

Chemical  Oxygen  Demand 

Dissoivod  Oxygen 

ParticulatB  Biogenic  Silica 

AvalableSifica 

Table  4-2 

Sediment  Medel  State  Variables  and  Fluxes 

State  Variabla 

SedimoM-Wster  FIor 

Temperature 

Particulate  Organic  Carbon 

Sediment  Oxygen  Demand 

SuMde/Methane 

Release  of  Chemical  Oxygen  Demand 

Particulate  Organic  Nitrogen 

Ammonium 

Ammonium  Flux 

Nitrate 

Nitiate  Fhtx 

Particulate  Organic  Phosphorus 

Phosphate 

Phosphate  Fhtx 

Particulate  Biogenic  SKca 

Avalabie  SiBca 

Sica  Fhtx 
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Mollifications  to  Water-Quality  Modal 

Witer>Celuiiiii  State  Variablas  and  Praeaaaas 

Initial  modifications  to  the  water-quality  model  will  be  die  addition  of 
particulate  inorganic  phosplmrus  and  zooplankton  to  the  suite  of  state 
variables.  Modifications  wUi  also  be  made  to  the  model  treatment  of 
dissolved  organic  conqKMinds  and  su^ioided  sediment. 

Particulate  Inorganic  Phoaidionis.  At  presmit,  litde  is  known  retarding  die 
nature  and  availability  of  particulate  inorganic  pbo^horus  in  the  bay  and 
tributaries.  The  pho^orus  may  be  tigfady  bound  in  mineral  complexes  and 
largely  unavailable  to  die  water  column  or  the  phosphorus  may  fr^y 
exchange  between  fractions  dissolved  and  loosely  sorbed  to  inorganic 
particles.  The  model  formulation  will  allow  for  bodi  exchangeable  and 
unexchangeable  forms  of  particulate  inorganic  phoqihorus.  Initial  guidance 
r^atding  these  forms  will  be  obtained  from  existing  data.  Furdier 
determinations  will  be  conducted  after  data  collected  as  part  of  die  tributary 
refinemaits  becomes  available. 

Zooptenkton.  One  or  two  zooplankton  groups  will  be  added  to  die  model. 
Detomination  of  die  number  of  groups  will  be  made  following  review  of 
available  data  and  investigation  of  computational  requirements. 

Dissolved  Organic  Compounds.  At  present,  the  model  rq[>resoits  labile  and 
refractory  fractions  of  particulate  organic  carbon,  nitrogen  and  phosphorus. 
Dissolved  organic  substances  are  not  divided  into  labile  and  refractory 
fractions,  however.  An  area  of  concern  in  die  next  phase  of  die  study  is  the 
potential  refractory  nature  of  nutrients  entering  the  bay  from  die  ocean. 

Present  rqiresentation  of  dissolved  organic  matter  will  be  revised  to  include 
labile  and  refractory  fractions.  The  revision  will  allow  inqiioved  investigation 
of  the  influence  of  oceanic  boundary  conditions. 

Suspended  Sediment.  Suspended  sediment  will  be  rqiresented  in  die  water- 
quality  model  to  account  for  interactions  with  particulate  inorganic  phosphorus 
aind  to  account  for  the  effect  on  light  extinction.  A  suspended-sediment  state 
variidile  is  currmitly  in  the  model.  Suqiended  sediment  is  treated  as  a 
conservative  substance  that  settles  at  a  specified  velocity.  At  present,  die 
suqimided-sediment  state  variable  is  employed  to  rqiresent  particulate  iron  and 
manganese  that  precipitate  following  autumn  turnover  of  anoxic  b<mom  water. 

The  mechanism  for  rqiresenting  suqiended  sediment  during  tributary 
refinemosts  will  be  determined  as  part  of  the  study.  Options  indude 
prediction  by  the  CH3D-WES  hyd^ynamic  model,  specification  based  on 
observations,  and  prediction  within  die  water-quality  model.  To  provide  for 
the  first  two  options,  the  water-quality  modd  code  wUl  be  modified  to  allow 
suspended  sediment  concentration  and  transport  to  be  read  in,  from  external 
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sources,  or  predicted,  as  in  the  present  code. 

Extension  of  Calibration.  Initial  nK}diftcations  to  the  water-quality  model 
will  be  performed  prior  to  availability  of  hydrodynamics  on  refined  tributary 
grids.  Preliminary  calibration  of  the  noodified  niodel  will  be  conducted  on  die 
existing  grid  using  existing  hydrodynamics  and  loads.  Transition  to  die 
refined  grids  will  be  completed  as  soon  as  the  revised  hydrodynamics  become 
available. 

The  Littoral  Zone 

The  littoral  zone  will  be  treated  as  a  distinct  region  in  the  wtder-quality 
model.  To  accommodate  the  littoral  zone,  the  model  code  must  be  c^g^  to 
allow  one  cell  thickness  along  the  shoreline  and  multiple  cell  dqpths 
elsewhere.  No  additional  changes  are  required  to  account  for  transport  in  the 
littoral  zone.  Cell  volume  and  transport  across  cell  faces  will  come  from  the 
hydrodynamic  model,  modified  to  include  die  littoral  zone  as  storage  area. 

Water  column  kinetics  and  state  variables  in  the  littoral  zone  will  be 
identical  to  elsewhere.  The  major  distinction  in  the  littoral  zone  will  be 
addition  of  a  submerged  aquatic  vegetation  (SAV)  conqionent.  Consensus, 
expressed  in  the  workshops,  is  that  SAV  cannot  exist  at  dqidis  greater  than 
2m,  even  under  ideal  circumstances.  Definition  of  the  littoral  zone  as  water 
inside  the  2m  contour  means  that  SAV  need  only  be  considered  within  the 
littoral  zone.  SAV  will  be  modelled  as  a  s^arate  submodel  which  will 
interact  widi  the  models  of  the  water  column  and  sediments.  Initial 
development  of  the  SAV  model  will  be  in  standalone  mode  widi  conditions  in 
the  water  column  and  sediments  specified  as  boundary  conditions. 
Development  in  standalone  model  will  focus  on  SAV  kinetics  without  the 
complications  and  resource  requirements  of  the  coupled  water-quality  and 
sediment  models.  A  pilot  for  this  process  was  development  of  the  sediment 
diagenesis  model  which  operates  in  standalone  or  coupled  mode. 

Formulation  of  the  SAV  model  will  be  determined  following  review  of 
SAV  models  currently  existing  for  die  bay  and  elsewhere.  Following  the 
recommendations  of  the  workshops,  initi^  formulation  of  the  SAV  model  will 
include  as  state  variables:  roots  and  rhizomes,  leaves,  epifauna,  and  grazers 
on  epifauna. 

Downstream  Boundary  Condition 

Nature  and  location  of  the  downstream  boundary  condition  in  the  water- 
quality  model  depends  upon  the  outcome  of  the  tests  conducted  with  the 
hydrodynamic  model.  If  mctension  of  the  grid  beyond  the  bay  mouth  is 
hydrodynamically  feasible,  the  downstream  boundary  of  die  water-quality 
model  will  be  at  the  limit  of  the  grid.  Concentration  boundary  conditions 
along  edges  of  the  grid  will  be  specified  based  on  observations  and  best 
estimates.  If  extension  of  the  grid  is  not  feasible,  other  altemadves  will  be 
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for  specificatioii  of  the  water-qiutlhy  boundary  condition.  One 
possibility  is  addition  of  a  few  well-mix^  cells  outside  the  bay  mouth. 

Nitrogen  and  phosphorus  concentrations  at  the  bay  mouth  have  decreased 
since  collection  of  the  1984*1986  obsovations  currently  onployed  in  the 
water-quality  nradel.  The  origin  and  nature  of  the  decrease  are  uncertain. 
Early  in  the  study,  sensitivity  to  the  recent  conditions  will  be  examined  in  the 
existing  model.  The  Monitoring  Subcommittee  will  be  requested  to  determine 
if  the  apparent  decreases  originate  in  laboratory  methodology. 

Modifications  to  Sediment  Model 

Tributary  Refinements 

Formulation  and  performance  of  the  sediment  model  in  tidal  freshwater 
will  be  enhanced  as  part  of  the  tributary  refinements.  In  the  present 
formulation,  sediment  model  parameters  are  arbitrarily  differentiated  into 
freshwater  and  saltwater  values  at  the  1  ppt  isohaline.  No  other  spatial 
differentiation  in  parameter  evaluation  is  possible.  No  obs^ations  in  tidal 
freshwater  and  few  observations  in  Virginia  tributaries  are  rq>resented  in  the 
present  calibration  data  base. 

The  minimum  suite  of  revisions  to  the  sediment  model  will  allow  for 
spatially-variable  specification  of  parameter  values.  An  attempt  will  be  made 
to  smoothly  transition  sediment  processes  from  saltwater  to  freshwater  instead 
of  differentiating  kinetics  at  an  arbitrary  isohaline.  Particulate  inorganic 
phosphorus  that  settles  from  the  water  column  will  be  routed  into  appropriate 
components  in  the  sediments.  The  model  will  be  calibrated  to  observations 
collected  in  appropriate  tributaries. 

At  present,  the  model  of  phosphorus  cycling  in  the  sediments  is  primitive. 
Phosphate  produced  by  diagenesis  sorbs  to  particles  or  is  released  to  the  water 
column.  More  detailed  modeling  of  phosphorus  cycling  b  desirable, 
especially  in  freshwater.  Improved  modeling  requires  addition  of  iron  and 
manganese  to  the  suite  of  sediment-model  state  variables.  Iron  and  manganese 
addition  will  also  improve  prediction  of  sediment  oxygen  demand  in 
freshwater.  The  feasibility  of  iron  and  manganese  modeling  will  be  examined 
and  the  two  metals  will  be  added  to  the  model  if  possible. 

Living  Resources 

Modeling  of  living  resources  also  necessitates  revisions  to  die  sediment 
model.  The  workshops  called  for  quantitative,  mechanistic  modeling  of 
benthic  macroinvertebrates.  Two  groups  of  macroinvertdirates,  suspmsion 
(heads-up)  and  d^sit  Oteads-down)  feeders,  will  be  added  to  die  suite  of 
sediment-model  state  variables.  Inclusion  of  macroinvertdirates  will  provide 
direct  quantification  of  the  biomass  of  this  important  food  source. 
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Rqiresentatioii  of  predation  on  algae  and  of  bioturbation  in  the  sediments  will 
also  be  enhanced. 

Other  modifications  to  the  sediment  model  will  be  addition  of  benthic 
algae  and  dissolved  organic  carbon  to  the  suite  of  state  variables  and 
acconunodation  of  nutrient  uptake  by  rooted,  submerged  aquatic  vegetation. 
These  features  are  required  to  model  the  littoral  zone. 

Sequence  to  be  Modelled 

Yoifc  River 

The  first  application  of  the  completely  revised  model,  including  the  littoral 
zone,  will  be  to  the  York  River.  Initial  (Vibration  will  be  on  a  sub-grid 
derived  from  the  grid  created  for  hydrodynamic  calibration  of  the  York  River. 
The  sub-grid  will  include  only  the  York  River  and  adjacent  bay. 
Hydrodynamic  boundary  conditions  at  the  open  edges  of  the  grid  will  be 
derived  from  the  hydrodynamic  model.  Water-quality  boundary  conditions 
will  be  derived  from  observations.  This  procedure  will  allow  a  large  number 
of  runs  to  be  made  rapidly.  The  York  River  will  be  calibrated  to  1985-1987 
data  in  one  continuous  run.  Following  calibration,  the  model  will  be  run  on 
the  complete  system-wide  grid  to  evaluate  performance  throughout  the  bay. 

Verification  of  the  York  River  model  will  be  conducted  employing 
observations  collected  In  1994.  The  verification  will  be  an  individual,  one- 
year  run.  Model  conditions  at  the  end  of  1987  will  be  employed  as  initial 
conditions  for  the  1994  run.  If  this  procedure  is  unsatisfactory,  consideration 
will  be  give  to  a  simulation  of  1985-1994.  Determination  of  Ae  grid  extent 
during  verification  (York  River  sub-grid  or  system-wide  grid)  will  be  made 
following  evaluation  of  model  performance  and  computational  requirements. 

James  and  Rappahannock  Rivers 

Calibration  of  the  James  and  Rappahannock  Rivers  will  commence  shortly 
after  calibration  of  the  York.  Calibration  activities  for  all  three  tributaries 
will  overliq).  As  with  the  York,  calibration  will  be  on  sub-grids  for  the  years 
1985-1987.  Verification  of  each  river  will  be  a  one-year  run  for  1994. 

Baltimore  Harbor 

Initial  calibration  of  the  Baltimore  Harbor  model  will  be  to  1985-1987 
observations.  The  monitoring  data  is  limited  to  one  station  in  the  Patiq[>sco 
and  one  in  Back  River.  Further  calibr^ion  and  verification  will  be  co^ucted 
against  the  1994-1995  data  as  it  becomes  available. 

Application  to  1994-1995  requires  loads  generated  by  a  watershed  modd 
simulation  of  those  years.  The  hydrodynamic  model  must  be  run  on  the 
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system-wide  grid  to  set  boundary  conditions  for  the  indq>endent  model. 
Employment  of  the  1994-1995  observations  within  the  time  frame  of  this 
project  depends  upon  the  availability  of  watershed  model  output  and  of  runoff, 
meteorological,  a^  other  data  needed  to  force  the  hydrodynamic  model.  We 
recommend  that  st^s  be  taken  now  to  ensure  required  inputs  are  available. 
Delay  in  delivery  of  required  inputs  will  result  in  increased  costs  and  delayed 
completion  of  the  project. 

Visualization 

The  present  model  produces  overwhelming  quantity  of  output.  Additional 
state  variables  and  grid  revisions  will  multiply  the  existing  output  volume.  A 
postprocessing  and  visualization  package  is  required  and  will  be  developed  as 
part  of  model  refinements.  Early  in  the  project,  incorporation  of  the  n^el 
into  the  Corps  FASTABS  system  will  be  explored. 

Scenarios 

Twenty  final  scenarios,  to  be  specified  by  the  Modeling  Subcommittee, 
will  be  run.  During  the  Three-Dimensional  Model  study,  decade-long  runs 
were  required  for  sediments  to  achieve  steady  state.  More  lengthy  runs  may 
be  required  in  the  tributaries  based  on  tributary  residence  time  and  local 
sediment  burial  rate. 

Procedures  for  running  the  scenarios  remain  to  be  determined.  One 
option  for  scenarios  in  Virginia  tributaries  is  to  run  scenarios  in  each  river 
using  the  system-wide  hydrodynamic  grid  with  emphasis  in  the  subject 
tributary.  A  second  option  is  to  combine  the  fine  grids  in  the  Virginia 
tributaries  into  a  single  system-wide  grid  with  fine  resolution  in  all  tributaries. 
The  first  option  rqiresents  system-wide  effects  with  minimum  computational 
requirements.  The  second  option  increases  computational  requirements  but 
produces  optimum  resolution  of  interactions  between  tributaries.  A  third,  fall¬ 
back  option  is  to  run  scenarios  on  tributary  sub-grids  using  boundary 
conditions  specified  from  the  current  three-dimensional  bay  model.  The 
scenario  procedure  will  be  determined  following  tests  of  mass-transport 
interactions  between  tributaries  and  evaluation  of  resources. 

Baltimore  Harbor  scenarios  will  be  run  on  the  fine  grid  using  results  of 
system-wide  scenarios  to  specify  boundary  conditions.  If  system-wide 
scenarios  are  not  run,  boundary  conditions  will  be  specified  by  an  alternate 
method,  perhaps  employing  the  existing  model. 

The  computational  requirements  to  run  decade-long  simulations  on  the 
system-wide  grid  with  all  model  modifications  are  immense.  Exact  time 
estimates  cannot  be  made  without  benchmarking  the  revised  model  but  order- 
of-magnitude  increases  in  computation  requirements  over  existing  scenarios 
are  feasible.  Scenario  execution  time  may  easily  increase  from  tiie  presrat  30 
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hours  CPU  time  on  a  Cray  Y>MP  to  300  hours.  A  combination  of  high-tech 
and  iowotedi  alternatives  to  connimption  of  6000  C7U  hours  (2S0  days)  for 
twenty  scenarios  need  to  be  investigated. 

Nalnr  Cnmgnlnn 

We  currently  have  available  a  Cray  C-90  udiich  executes  finty-percent 
foster  than  the  Y-MP.  Three  years  from  now,  still  foster  conqniters  will  be 
available.  Alternate  conqmters  may  also  be  comidered.  Early  in  the  project, 
investigation  of  current  model  performance  on  a  massive  parallel  processor 
will  be  conducted. 

Revised  Code 

One  q[)tion  to  reduce  computation  time  is  to  update  kinetics  and  sediment 
fluxes  at  less  frequent  intervals  than  eadi  time  st^.  No  need  exists,  for 
example,  to  recompute  SOD  every  hour.  Daily  updates  will  likdy  suffice. 

The  benefits  of  less-frequent  updates  will  be  investigated.  A  ”fix*  to  speed  up 
die  sediment  model  time  to  su^y  state  was  proposed  during  the  Three- 
Dimensional  Model  Study  and  may  be  inqilanented  in  tributary  refinement 
scenarios. 

Revisions  described  above  can  be  impiemented  by  the  WQM  team.  WES 
wUl  also  solicit  advice  and  services  of  experts  in  optimizing  die  code. 
Optimization  will  be  conducted  in  two  phases.  One  round  of  optimization  will 
be  conducted  at  commencement  of  die  project.  The  optimized  code  will  be 
nearly  identical  to  the  existing  code.  Revisions  wOl  follow  die  pattern  set  by 
consultants  to  the  EPA  vdio  halved  the  execution  time  of  the  present  model. 

A  second  phase  of  qitimization  will  be  conducted,  prior  to  pr^uction  scenario 
runs.  At  this  phase,  the  qitimized  code  will  be  close  to  die  final  code 
developed  for  tributary  reflnemoits. 

Reduced  Grid 

The  system-wide  grid  can  be  ov^aid  to  reduce  die  nundier  of  cells  away 
from  r^ions  of  interest.  Anodier  alternative  is  to  run  scenarios  on  die  James 
and  Rqipahannock  enqiloying  die  subgrids  developed  for  rqiid  tributary 
calibration.  Boundary  conditions  would  be  specified  by  runs  on  only  one  of 
the  system-wide  grids,  e.g.  die  York  River  grid.  Yet  another  alternative  is  to 
run  scenarios  on  a  grid  whidi  includes  the  lower  half  of  die  bay  only. 
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Suspended-Sediment  Activities 


Suspended-Sediment  Budget 

The  first  stq>  in  modding  subtended  sediments  will  be  investigation  of  die 
suspended-sedin^  budget  of  die  system.  Hie  primary  objective  is  to 
determine  if  sources  of  sediment  are  known  sofficiendy  to  allow  for  a  more 
advanced,  mechanistic  sediment  model.  Additional  objectives  are  to  indicate 
the  fraction  of  the  sediment  inputs  diat  are  managedile  and  to  investigate  die 
feasibility  of  employing  die  water-quality  model  as  a  suqiended-sediment 
model. 

The  existing  wata--quality  model  will  be  enqiloyed  in  a  mass-balance 
accounting.  The  accounting  will  indicate  if  obs^^  suqiended  sediment 
concoitrations  can  be  computed  based  on  known  sediment  loads.  Fall-line 
loads  will  be  alternately  d^ved  from  n^resskm  on  observations  and  from  die 
watershed  model.  Below-fidMine  noiqxiint-source  loads  will  be  derived  fipom 
the  watershed  model,  if  av^able.  Other  bdow-fall-line  sources,  notably  bank 
erosion,  will  be  derived  from  best  available  information.  Oceanic  boundaries 
will  be  obtained  from  observations.  Net  settling  will  be  specified  based  on 
measures  of  long-term  burial. 

Two  classes  of  sediment  will  be  considered,  organic  and  inorganic. 
Observed  organic  solids  will  be  derived  from  (rt»ervations  of  particulate 
organic  carbon  in  die  data  bue.  Observed  inorganic  solUs  will  be  determined 
as  observed  total  su^iended  solids  less  organic  solids.  Predicted  organic 
solids  will  be  detemtined  from  presentiy  modeled  particulate  organic  carbon. 
Predicted  inorganic  solids  will  be  confuted  by  mass  balance  using  die  existing 
modd  suspended  solids  state  variable. 

Concentrations  and  spatial  distributions  of  predicted  and  observed  solids 
will  be  compared  in  the  mainstem  and  tributaries.  Predicted  concentrations 
will  indicate  our  ability  to  "balance  die  budget".  Spatial  distributic  ^’ill 
provide  a  first  look  at  die  feasibility  of  predicting  pbenomena  sudi  c  die 
turbidity  maximum  using  die  water-quality  siodd. 


Chasers  SOTptnSad Sadliiitnt AciMSh 


5-1 


Wind  Resuspension  of  Sediments 

Investigations  indicate  Aat  submerged  aquatic  vegetation  (SAV)  is 
sensitive  to  short-term  fluctuations  in  light  extinction  induced  by  fluctuations 
in  suspended  sediment.  A  prinury  agent  of  sediment  fluctuations  is  wind 
resuspension.  Consequently,  a  medianistic  rq>resentation  of  wind 
resuspension  of  sedimoits  will  be  incorporated  into  die  modd  of  the  littoral 
zone.  The  resuspension  modd  will  be  based  on  existing  data,  primarily 
collected  by  VIMS  in  the  York  River.  The  relationship  of  light  extinction  to 
suspended  solids  will  be  based  on  existing  data  and  on  observations  collected 
during  enhanced  tributary  monitoring.  The  sufficiency  of  these  data  bases  for 
modd  cdibration  and  vmfication  is  unknown.  Additional  studies  to  measure 
wind  resuspension  in  Virginia  tributaries  and  Baltimore  Harbor  and  to 
examine  die  rdationship  of  light  and  solids  are  under  consideration.  Data 
devdoped  in  these  studies  will  be  enqiloyed  as  it  becomes  available. 


5-2 


CliaplwS  Sinpanditf -Stdimwit  AmMSm 


Chapter  6 

Coordinating  Efforts 


Devdopmeot  and  calibration  of  the  refined  tributary  nd  livim-reaoarce 
models  will  be  nqiported  by  data  in  the  data  base  maintained  by  the  CBPO, 
by  ou^Mit  fnmi  die  Watershed  Modd,  and  by  observations  conducted 
qiedficdly  for  die  projed.  The  present  duqiter  sunmuuizes  the  coordfoadag 
efSMts  required  for  successful  completion  of  the  study. 

Existing  Data 

Upon  commencement  of  die  study,  a  detailed  request  for  data  wfll  be 
{wesemed  to  the  CBPO.  Most  signified  information  to  be  requested  is  listed 
bdow. 

Monitoring  Data 

Observations  collected  in  die  mdnstan,  trilHitaries,  and  fdl  lines  by  the 
monitoring  program  will  be  siqiplied  by  die  CBPO  to  WES  in  the  form  of 
SAS  datasets.  In  die  event  SAS  datasets  cannot  be  supplied,  anodier  format 
compatible  widi  WES  cennputers  wfll  be  provided.  Data  must  be  supplied  for 
1985-1987  and  1994-1$195,  at  least.  Provision  of  all  monitoring  data,  1984  to 
date,  is  preferable,  however.  The  data  wfll  be  avaflaUe  for  potential  modd 
extension  and  for  purposes  odier  dum  modd  data  comparison. 

Zooplankton  and  Maerobonthos 

Zoc^lankton  and  bendiic  macrobendios  ue  monitored  by  die  CBPO  but 
are  not  routindy  provided  as  part  of  die  monitoring  data  base.  Observatkms 
wfll  be  provided  to  WES  in  SAS  data  sets  or  odier  conqiadble  format.  To  be 
of  use,  observations  must  be  in  units  of  biomass,  prefenbly  as  carbon.  If 
observatkms  are  suiqilied  as  number  of  organisms,  suitable  conversion  foctors 
to  biomass  must  be  siqiplied  as  wdl. 


Siibni«iv«d  Aquatic  Vagatatioa 

Submerged  aquatic  v^etation  acreage  and  density  wUl  be  provided  ftv 
years  of  record.  Information  will  be  provided  for  ea^  modd  cdl  based  on 
cdl  coordinates  supplied  by  WES  to  foe  CBPO.  Since  several  grids  will  be 
enq>loyed,  multiple  data  requests  will  be  completed. 

Watershed  Model 

Water-Qualhy  Calibration  Yean 

Application  of  foe  water-quality  modd  to  foe  years  198S-1S187  requires 
execution  of  foe  watershed  m^d  for  foe  same  years.  Prior  linkage  of  foe 
watershed  modd  to  foe  water-quality  modd  employed  loads  specified  on  a  bi¬ 
weekly  or  monthly  basis.  Bdow-foll-line  loads  were  distribute  uniformly 
along  the  length  of  each  tributary.  Both  foe  Three-Dimensional  Modd 
calibration  rq)ort  (Cerco  and  Cole  1994)  and  foe  first  workshop  (HydroQual 
1993)  called  for  increased  temporal  and  spatial  resolution  of  loads  to  foe 
water-quality  modd.  Loads  can  be  readily  iiqmt  to  foe  modd  on  a  daily 
basis.  The  worlq)lan  for  Phase  in  Waterfoed  Modd  efforts  calls  for  refined 
spatial  detaU  but  not  in  Virginia  tributaries.  As  an  alternative,  loads  in  bdow- 
^1-line  watersheds  will  be  dlocated  into  major  sub-basins.  Loads  outside 
major  sub-basins  will  be  uniformly  distributed  dong  foe  lengfo  of  eadi 
tributary.  Allocation  of  load  into  sub-basins  will  be  conducted  by  foe  CBPO. 
WES  will  assign  sub-basin  loads  to  appropriate  model  cdls. 

Water-Quality  Verification  Years 

Application  of  foe  water-qudity  modd  to  observations  collected  in  1994- 
199S  requires  execution  of  foe  watershed  modd  for  foe  same  years.  In  view 
of  foe  time  required  to  assemble  watershed  modd  iiq>ut  data,  plans  should  be 
made  immedu^y  to  acquire  foe  data.  Otherwise,  dternatives  to  aiq>loyment 
of  foe  watershed  modd  need  to  be  devdoped. 

The  1985-1987  loads  will  be  initially  produced  by  Phase  m  of  foe 
watershed  modd.  The  1994-1995  loads  wUl  be  produced  by  Phase  IV  of  foe 
watershed  modd.  Cdibrating  foe  water-qudity  modd  to  Phase  in  loads  and 
verifying  with  Ifoase  IV  loads  presents  potentid  problems.  Results  from  foe 
Ifoase  IV  version  of  foe  waterfoed  modd  should  be  carefully  diecked  by  foe 
CBPO  for  consistency  with  Phase  DI.  If  substantid  differences  exist,  1985- 
1987  loads  must  be  produced  from  Phase  IV  and  iiq>ut  to  foe  water-qudity 
modd. 

Point-Source  Characterization 

As  part  of  foe  effort  to  refine  loading  information  to  foe  modd,  effluent 
characteristics  will  be  measured  at  42  facilities  in  Virginia  and  12  in 
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Loads  from  each  point  source  will  be  ocmqNited  by  the  CBPO  on  a 
mmitfily  basis  fiw  die  raodd  spplicatkm  period,  bfoadily  loads  will  be 
siqiplied  to  WES  as  SAS  dalasi  or  in  oonqiatible  format. 

Supporting  Observations 

ObservatkMis  beycmd  the  routine  monitorii^  program  will  be  ocmducted  in 
the  tributaries  and,  to  a  lesser  extern,  in  the  mainsrem  and  Pottmiac. 

Cdlectkm  in  Virginia  tributaries,  die  mairndem,  and  Potomac  wfll  be 
conducted  primarily  in  1994.  Cdlection  in  Baltimore  Harbor  wfll  be  in  1994 
and  1995. 

Viifinfai  Tributaries 

Addidonal  Stathms.  Monitoring  stations  wfll  be  added  along  die  longitudinal 
and  lateral  axes  of  die  tributaries.  Additional  lateral  stations  will  be  located 
widiin  die  littoral  xone,  halfway  between  the  2  m  contour  and  die  shordine. 

Add thmal  Substances,  Particulate  inorganic  jflioqdiorus  and  pardcuiate 
biogenic  silica  wfli  be  added  to  die  suite  of  nxmitoring  variables.  These  wfll 
be  analyzed  in  dl  frdMine  saiqiles.  The  additkmd  substances  wfll  be 
andyzed  at  dl  fahstream  statkms  but  for  only  hdf  the  surveys. 

SOME  Measures.  Measures  of  sedimem  oxygen  demand  and  nutrient  flux 
wfll  be  conducted  in  two  of  the  tributaries,  the  James  and  York.  Measures  in 
tidd  freshwater  wfll  be  conducted  in  the  James  only.  Measures  wfll  also  be 
conducted  that  cmnpate  fluxes  in  die  channd  versus  littord  zone.  Littord 
zmie  measures  wfli  be  cmiducted  in  li^t  and  dark  to  evduate  effects  of 
bendiicdgae. 

Bahhnora  Harbar 

Adtfdonal  Stattons.  The  nnndier  of  monitoring  stations  in  die  Pd^isoo 
River  wfll  be  increased  from  one  to  nineteen.  Mtmitoring  stations  in  the  Back 


River  will  be  increased  from  one  to  four. 

Additional  Substances.  Particulate  inorganic  pho^danus  and  particulate 
biogenic  silica  will  be  analyzed  at  all  stations.  Sanqiles  will  be  collected  in 
only  a  portion  of  die  surveys,  however. 

SOME  Measures.  Measures  of  sedimmit  oxygen  demand  and  nutrient  flux 
will  be  conducted  at  diree  stations  in  Baltimore  Harbor. 

Sediment  Resuspension.  A  program  is  prt^xised  to  collea  sediment  and 
hydrographic  data  required  to  devdop  and  Vibrate  a  fully^iredictive 
suspended-sediment  tran^rt  model. 

Mainstoro  Bay  and  Potomac 

Additional  Substances.  Particulate  inorganic  phosidiorus  and  particulate 
biogenic  silica  will  be  analyzed  at  all  mainstem  stations  and  de^hs  at  wfaidi 
nutrients  are  analyzed  in  the  monitoring  program.  Particulate  inorganic 
phosphorus  and  particulate  biogenic  silica  will  be  analyzed  at  three  stations  in 
the  Potomac.  In  both  systems,  sanqiles  will  be  collected  in  only  a  portion  of 
the  surveys. 

Potential  Additional  Observation  Programs 

The  workshops  and  subsequent  nmetings  indicated  numerous  data 
collection  efforts  of  value  to  the  tributary  refinements.  Resources  available 
when  the  program  of  extended  tributary  monitoring  was  planned  did  not  allow 
for  all  data  requirmnoits  to  be  fulfilled.  Odier  data  needs  were  not  obvious 
vdien  the  exte^ed  monitoring  was  planned.  A  list  of  recommended  programs 
has  been  prq>ared  by  the  CBPO  and  will  be  considered  by  die  Monitoring 
Committee.  These  include: 

Detailed  study  of  the  relation  of  light  absorbance  and  scattering  to 
dissolved  and  suspended  solids. 

Measurement  of  sedirnem  resuspension  by  wind  in  existing  and 
potential  SAV  locations. 

Additional  SOME  measures  in  tidal  freshwater. 

Measure  of  labile  and  refractory  portions  of  dissolve  organic  matter  at 
the  ocean  boundary. 

Collection  of  continuous  salinity  observations  for  validation  of 
exdiange  between  littoral  zone  and  main  channel. 
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Living-Resources  Woitshops 

The  Living  Resources  Subcommittee  of  die  CSiestgieake  Bay  Program  has 
sponsored  numerous  data  collection  and  ecosystem  moving  efforts. 
Observations  and  knowledge  developed  during  these  efforts  are  invaluable  to 
die  living-resource  modeling  proposed  for  the  present  study.  Annual 
workshops  will  be  conducted  to  promote  cooperation  and  information 
exchange  betwemi  investigators  qionsored  by  the  Living  Resources 
Subcommittee  and  participants  in  the  presnt  study.  The  woriohops  will  be 
organized  by  die  CBPO  or  one  of  its  contractors. 
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